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ABSTRACT 

We present the result of our low-luminosity quasar survey in the rcdshift range of 4.5 < z < 5.5 in 
the COSMOS field. Using the COSMOS photometric catalog, we selected 15 quasar candidates with 
22 < i! < 24 at z ~ 5, that are ~ 3 mag fainter than the SDSS quasars in the same redshift range. 
We obtained optical spectra for 14 of the 15 candidates using FOCAS on the Subaru Telescope and 
did not identify any low-luminosity type-1 quasars at z ~ 5 while a low-luminosity type-2 quasar at 
z ~ 5.07 was discovered. In order to constrain the faint end of the quasar luminosity function at 
z ~ 5, we calculated the la confidence upper limits of the space density of type-1 quasars. As a 
result, the la confidence upper limits on the quasar space density are <fr < 1.33 x 10~ 7 Mpc~ 3 mag -1 
for -24.52 < A/i 450 < -23.52 and $ < 2.88 x lCr f Mpc~ 3 mag" 1 for -23.52 < M U5a < -22.52. 
The inferred la confidence upper limits of the space density are then used to provide constrains on 
the faint-end slope and the break absolute magnitude of the quasar luminosity function at z ~ 5. 
We find that the quasar space density decreases gradually as a function of redshift at low luminosity 
(M1450 ~ —23), being similar to the trend found for quasars with high luminosity (M1450 < —26). 
This result is consistent with the so-called downsizing evolution of quasars seen at lower redshifts. 
Subject headings: cosmology: observations — quasars: general — surveys 



1. INTRODUCTION 

The evolution of supermassive black holes (SMBHs) is 
now regarded as one of the most important unresolved 
issues in the modern astronomy, after the discovery of 
the galaxy-SMBH "co-evolution" inferred from, e.g., a 
tight relationship between the mass of SMBHs and their 
host bulges (e.g.. iMarconi fc Hunt! 12001 IHaring fc Rixl 
12004 iGultekin et al.ll2009fl . Measuring the whole shape 
of the quasar luminosity function (QLF) is particularly 
important to understand how the SMBHs grew, since it 
is highly dependent on some key parameter s of SMBHs 
such as the growth timescale of SMBHs (e.g. JEnoki et al.l 
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12001 . 

The QLF at z < 3 has b een well quantified over a 
wide luminosity range (e.g., iCroom et al.l 120091) and is 
best r epresented by a double power law (e.g.. lBovle et al.l 
[19881: iPeil fl995h . Recently, the f aint end of the QLF 
at z ~ 4 has been measured (IGlikman et al.l 120101 : 
llkeda et alJl20Tfl IGlikman et al.l 120111 ) and is also best 
represented by a double power law. More interestingly, 
recent studies on the optical QLF show that the space 
density of low-luminosity active galactic nuclei (AGNs) 
peaks at a lower redshift than that of mor e luminous 
AGNs (jCroom et al.l l2009t llkeda et al.l 120111 ). This re- 
sult can be interpreted as AGN (or SMBH) downsizing 
evolution, since the brighter AGNs tend to harbor the 
more massive SMBHs if the dispersion of t he Eddington 
ratio of quasars is not very large (see, e.g., iTrump et al.l 
120111) . The A GN downsizing has been also reported by 
X-r av surveys (|Ueda et a l. 2003; Hasi nger et al.l2005t see 
also iBrusa et al.ll2009l and lCivano et al.ll2011l ). However , 
the physical origin of the AGN downsizing is totally un- 
clear, that makes h igh-z low-luminosity quasar surveys 
more important (see Fanida kis et al.l l2012. for theoretical 
works on the AGN downsizing evolution). 

Recently, some low-luminosity qu asar survey s have 

been perform ed at z > 5 (IC00I et al.l 120061 : 

iMahabal et aHl2005h . IC00I et al.1 (|2006|) identified three 
quasars at z > 5 with z' < 22 and included a quasar 
at z = 5.85 with z' = 20.68, in the AGES survey which 
covers 8.5 deg 2 . Jiang et al. (2008) also identified five 
new quasars at z > 5.8 with 20 < z 1 < 21 in the Sloan 
Digital Sky Survey (SDSS) deep stripe which covers 260 
deg 2 . The space density of quasars at z ~ 6 which is 
calculated by the result of Cool et al. (2006) is about six 
times larger than the result of Jiang et al. (2008). This 
large discrepancy may be caused by the small survey area 
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of Cool et al. (2006). iMahabal et all ((2005ft identified a 
very faint quasar at z = 5.70 with z' = 23.0 in the tota l 
quasar survey area of 2.5 deg 2 . M ahabal et all §005) 
mentioned that the surface density of quasars at similar 
redshifts is roughly consistent with previous extrapola- 
tions of the faint end of the QLF. In this way, some low- 
luminosity quasars have been discovered although the 
faint end of QLF is not determined exactly, due to the 
lack of low-luminosity quasars. 

At z > 6, a number of luminous qu asars have been 
found up to z ~ 6.5 by the SPSS (e .g., iFan et all [2001 
iGotol [2001 Uiang et al.l [200l [2009ft and the Canada- 



Franc e High-z Quasar Su r vey (CFHQS; Willott et al.l 
120071 IWillott et all [200l IWillott et all 12010ft . Re^ 
cently, a luminous qu asar at z = 7.085 has been found 
(jMortlock et aljfeoill ) through the data obtained by the 
United Kingdom Infrared Telescope In frared Deep Sky 
burvey qjKIDSS: lLawrence et al.ll2007h . Although some 
low-luminosity quasars at z > 5 have been discovered 
as mentioned above, the faint-end slope of the z > 5 
QLF is still very poorly constrained. Consequently it 
is not understood how low-luminosity quasars evolve at 
high redshifts, or if the AGN downsizing evolution is also 
seen in the earlier universe. Since the number density of 
low-luminosity quasars is expected to be much higher 
than that of high-luminosity quasars, the whole picture 
of SMBH evolution cannot be understood without firm 
measurements of the number density of low-luminosity 
quasars at such high redshifts. 

Motivated by these issues, we have searched for 
low-luminosity quas ars at z ~ 5 in the COSMOS field 
(|Scoville et al .112007ft . In Section 2, we describe the data 
and method that were used for the photometric selection 
of quasar candidates. In Section 3, we report the 
results of the follow-up spectroscopic observations. In 
Section 4, we describe how we estimate the photometric 
completeness to derive the QLF. In Section 5, we present 
the upper limits of the QLF at z ~ 5 and briefly discuss 
it. Throughout this paper we use a A cosmology with 
fi TO = 0.3, Ha — 0.7, and the Hubble constant of H = 
70 km s _1 Mpc" 1 . 

2. THE SAMPLE 
2.1. The Cosmic Evolution Survey 

The COSMOS is a treasury program on the Hubble 
Space Telescope (HST) and comprises 270 and 320 orbits 
allocated with HST Cycles 12 and 13, respectively (Scov- 
ille et al. 2007; Koekemoer et al. 2007). The COSMOS 
field covers an area of 1.4° x 1.4° square which corre- 
sponds to 2 deg 2 , centered at R.A. (J2000) = 10:00:28.6 
and Dec. (J2000) =+02:12:21.0. 

We use the official C OSMOS photometr ic redshift 
catalog for phot ometry (jllbert et al.l [2009J see also 
ICapak et al.ll2007l ) to select the quasar photometric can- 
didates at z ~ 5. This catalog covers an area of ~ 2 
deg 2 and contains several photometric measurements. 
Specifically in this paper, we use the it*-band 3 di- 
ameter aperture apparent magnitu de measured on the 
image obtained with MegaCam ([Boulade et al.l 12003ft 
on the Canada- France-Hawaii Telescope (CFHT), and 
the 3 diameter aperture apparent magnitudes of the 
g'-, r'-, i'-, and z'-bands (|Taniguchi et al.l [ 20 07) mea- 
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Fig. 1. — Two color diagram of r' — %' and i' — z', that is used for 
our quasar selection. Green plus symbols denote point sources with 
22 < t'(MAG_AUT0) < 24. Blue an d magenta sq uares denote colors 
of M-type stars and K-type stars (Pickles 1998), respectively. The 
red line is the median track of the model quasar colors. The red 
error bars show the standard deviation of the r' — i' and V — z 1 
colors in our model quasar spectra. The black solid line shows our 
photometric criteria used to select quasar candidates at z ~ 5. 

sured on the image obtaine d with the Subaru Suprime- 
Cam (Mivaza ki et al.l [ 20 02). and the z'-band total ap- 
parent magnitude (MAG_AUTD measured by SExtractor; 
iBertin fc Arnoutsl Il996[ ) whose measurement is also 
based on the Suprime-Cam i'-band image. 

The 5cr limiting AB apparent magnitudes are u* = 
26.5, g' = 26.5, r> = 26.6, i' = 26.1, and z' = 25.1 (3" di- 
ameter aperture). Since we also use the Advanced Cam- 
era for Surveys (ACS) catalog (|Koekemoer et al.l 120071 : 
iLeauthaud et al.l 120071 ) to separate galaxies from point 
sources, our survey area is restricted to the area mapped 
with ACS on HST (1.64 deg 2 ). Note that all of the 
data in the official COSMOS photometric redshift cat- 
alog overlaps the entire ACS field. 

2.2. Quasar Candidate Selection 

Quasars at z ~ 5 show the Lyman break in their spec- 
tral energy distribution (SED) that falls between the 
wavelengths of the r' and i' filters, making their 



color redder than their 



color. We utilize this 



property to select candidates of low-luminosity quasars 
at z ~ 5. Here typical quasar colors as a function of 
redshift are necessary to define reliable color-selection 
criteria for quasars at z ~ 5. Therefore we generate 
model quasar spectra following t he procedure generally 
adopted in previous studies (e.g.. lFanl|1999b IHunt et"aLl 
12001 IRichards et al.H200l ISiana et al.ll2008l ). and derive 
the g'—r', r'—i', and i'—z' colors of the model quasars at 
redshifts from to 6. In this procedure, we adopt the typ- 
ical power-law slope (a u = 0.46, where f u oc i>~ au ), Lya 
rest-frame equivalent wi dth (EWp = 90 A), and t ypical 
emission- line flux ratios (fVanden Berk et alj |200l!). The 
effects of the intergalactic absorption by the neutral hy- 
dr ogen we r e corr ected by adopting the extinction model 
of iMadaul (|1995ft . Our simulated colors of the model 
quasars are shown in the r' — i' versus i' — z' diagram 
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TABLE 1 

Properties of the quasar candidates at z ~ 5 



Number 


R.A. 


Decl. 


i' (MAG.AUT0) 


r> - V 


V - z' 


Exp.Time a 




(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


(sec) 


1 


150.69131 


1.637161 


23.40 


2.34 


0.67 


2400 


2 


150.45275 


1.669653 


23.48 


2.04 


2.69 


3000 


3 


150.17448 


1.629074 


23.76 


1.67 


2.11 


1800 


4 b 


150.64917 


1.816186 


23.39 


3.44 


0.82 


- 


5 


149.87082 


1.882791 


23.98 


1.26 


0.17 


2400 


6 


149.85403 


1.823611 


23.97 


2.58 


0.87 


2400 


7 


149.78245 


2.221621 


23.96 


3.19 


1.13 


1800 


8 


149.69804 


2.283260 


23.67 


2.04 


0.44 


1800 


9 


150.56861 


2.317432 


23.98 


4.09 


1.26 


1800 


10 


150.05481 


2.376726 


23.89 


1.09 


0.25 


2700 


ll c 


149.78381 


2.452135 


23.70 


1.35 


0.26 


7200 


12 


150.16401 


2.549605 


23.31 


1.96 


0.61 


2400 


13 


149.96443 


2.473646 


23.93 


1.21 


0.21 


2700 


14 


150.66035 


2.786445 


23.51 


1.93 


0.57 


2400 


15 


149.59161 


2.659749 


23.16 


2.26 


0.77 


2400 



"Total on-source exposure time in the FOCAS spectroscopic observation. 
b This quasar candidate was not observed. 
c Type-2 quasar at z = 5.07. 



(Fig. 1). Note that the dispersions in the power-law 
slope and EWs of quasars are also taken into account 
when we calculate the photometric completeness of our 
survey (see Section 4). 

We select our quasar photometric candidates at z ~ 5 
using the following criteria: 

22 < i'(MAG_AUT0) < 24, (1) 
i' - z' < 0.45(r' - i') - 0.24, (2) 

u* > 27.49, (3) 

,g'-r'>1.3, (4) 

and 

r'-i'> 1.0, (5) 

The criterion (2) is used to select quasars efficiently with- 
out significant contamination from stars (especially from 
M0 to M6, see also Fig. 1), taking the color distribu- 
tion of stars and model quasars into account. To re- 
move possible foreground contaminations further, we in- 
troduce the additional criteria (3), (4), and (5). These 
latter two color thresholds arc adopted by taking empir- 
ical color distribution s of quasars at z ~ 5 into account 
(|Richards et al.|[200l . 

Here we comment on our point-source criterion based 
on the HST image (F814W, see Koekemoer et al. 2007), 
whose spatial resolution is 0'.'09 in FWHM (that cor- 
responds to ~ 0.6kpc at z = 5). Since the size 
of high-z quasar h ost galaxies is larger than 1 kpc 
(|Targett et al.lf20T2h . it is expected that the quasar host 
galaxies are spatially resolved in the ACS images. There- 
fore it would be possible some quasars could be ex- 
cluded from the sample of our quasar photometric can- 
didates. However, at z ~ 1, some earlier works show 
that the host galaxy of quasars with a similar abso- 
lute magnitude to our targe ts is typically ~ 2 mag 
fainter than their nucleus ( e.g., lMcLeod fc McLeodll2001l : 
iMcLeod fc Bechtoldl [20091 ). The brightness difference 



may be even larger at higher redshifts, because the typ- 
ical Eddington ratio of luminous quasar s is roughly con- 
stant at z ~ 1 — 4 (jTrump et al.ll2009T ) while the mass 
ratio of SMB Hs to host galax i es is higher at higher 
redshifts (e.g . iTreu et all 120061 IWoo et al.ll2fM [20061 : 
iMerloni et all 120101: iBennert et al.ll2010L 1201 ID . There- 
fore we conclude that quasars explored in this work 
should be recognized as point sources in the HST image. 
To distinguish the galaxies and point sources, Leauthaud 
et al. (2007) used the SExtractor parameter MU_MAX 
(peak surface brightness above the background level) and 
MAG_AUT0 (see Fig. 5 of Leauthaud et al. 2007). Because 
the fact that the light distribution of a point source varies 
with magnitude. Therefore we can distinguish the ex- 
tended objects and point sources by using the MU_MAX 
and the MAG_AUT0. 

Accordingly we removed 23 spatially-extended objects 
satisfying the criteria (l)-(5) based on the classification 
by Leauthaud et al. (2007). As a result, we obtain a 
sample of 15 quasar candidates among 7318 point sources 
with 22 < i'(MAG_AUT0) < 24. The selected candidates 
are listed in Table 1. 

3. SPECTROSCOPIC OBSERVATION 

The spectroscopic follow-up observations of the quasar 
candidates were carried out at the Subaru Telescope 
with the Faint Object Cam era and Spectrograph (FO- 
CAS; |Kaihi^^aiLil]|l002) on 7-11 January 2010. We 
used the 300 grating with the S058 filter, whose wave- 
length coverage is 5800A < A obs < lOOOOA. We used 
a O'.'8-width slit, resulting in a wavelength resolution of 
R ~ 700 (Av ~ 430 km s _1 ) as measured by night sky 
emission lines. The typical seeing size was ~ 0'.'7. Due 
to the limited observing time, we observed 14 of the 15 
candidates; the object No. 4 in Table 1 was not observed. 
The individual exposure time was 600 - 900 sec, and the 
total exposure time was 1800 - 7200 sec for each object. 

Standard data reduction procedures were performed 
by using IRAF. After the sky subtraction, we extracted 
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Fig. 2. — Reduced spectrum of object No. 11 in Table 1 (upper 
panel) and typical sky spectrum (lower panel). Dotted lines show 
the expected wavelengths of quasar emission lines: Lya A1216 and 
C iv A1549. 

one-dimensional spectra with an aperture size of l'/8 and 
the relative sensitivity calibration was performed using 
the spectral data of a spectrophotometric standard star, 
Feige 34. The spectra of 14 objects were then flux- 
calibrated using the i'-band photometric magnitude of 
these objects. We found that one spectrum shows only 
narrow Lya emission lines at A Q b s = 7381 A(z ~ 5.07 and 
A^fwhm ~ 800 km s _1 ) without any high- ionization 
lines such as C iv (Fig. 2). Since this object is de- 
tecte d in the X-ray ban d by the Chandra-COSMOS sur- 
vey (jElvis et all 120091 CID-2220) and its X-ray lumi- 
nosity is 3xl0 44 erg/s in the 2-10 keV rest frame band 
(jCivano et alJl201lU2012l ). it can be classified as an AGN. 
We can classify this object as a Type 2 AGN based on the 
upper limit available for the X-ray hardness ratio consis- 
tent with mild obscuration (NH<5xl0 23 cm -2 ) together 
with its Ly a spectral profile. Although the X-ray hard- 
ness ratio is not available for this object, we classify this 
object as a type-2 AGN based on its Lya spectral pro- 
file. Therefore we conclude that no type-1 quasars are 
identified in our spectroscopic follow-up campaign. The 
spectra of the remaining 13 objects are consistent with 
those of Galactic late- type stars, and an example of these 
spectra is shown in Fig. 3. 

4. COMPLETENESS 

Quasar surveys are generally not perfectly complete 
due to various factors such as photometric errors and 
intrinsic variations in the spectra. Therefore, to derive 
the QLF acculately, the survey completeness needs to be 
estimated as a function of the quasar redshift and ap- 
parent magnitude. We derive the completeness by mod- 
eling quasar spectra, in a similar way as described in 
Section 2.2. Here we also take into account the intrinsic 
variation in the continuum slope and EWs of the emis- 
sion lines. We assume a Gaussian distribution of the 
power-law slope a„ (f v oc v~ a ") and Lya EWs, with 
means of 0.46 and 90A (the same as those in Section 2.2), 
and a standard deviation of 0.30 and 20A, respectively 
(IVanden Berk et al.|200lHFrancidH99l^lHunt et al.ll2004J: 
Glik man et alJl2O10l T" We include emission lines whose 
flux is larger than 0.5% of the Lya flux, given in Table 
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Fig. 3. — Example of the late type star spectrum (upper panel; 
No. 2 in Table 1) and typical sky spectrum (lower panel). 

2 of Vanden Berk et al. (2001). The emission-line ratios 
are assumed to be the same for all model quasars (i.e., 
scaling to the Lya strength) . We also include the Balmcr 
continuum an d Fe I I features by using the template of 
iKawara et al.l (|1996l ). We create 1000 quasar spectra at 
each Az = 0.01 in the redshift range of < z < 6. The ef- 
fects of intergalactic absorption by neutral hydrogen were 
corrected by adopting the extinction model of Madau 
(1995). Then, we calculated the colors of the model 
quasars in the observed frame. We compared the colors 
of the simulated quasars with the empirical quasar colors 
(SDSS DR7) to check whether the simulated quasar col- 
ors are consistent with the empirical quasar colors. Fig. 
4 shows the comparison between the empirical and sim- 
ulated quasar colors of g* — r* , r* — i*, and i* — z* by 
using the transmission curves of the SDSS filters. Note 
that the dispersion of the simulated quasar colors is sys- 
tematically smaller than the dispersion of the observed 
quasar colors, through the average color is consistent be- 
tween the simulated and observed quasar colors. This is 
because the simulated colors presented here do not take 
photometric errors into account. The photometric errors 
are properly taken into account when the completeness 
is calculated, as described below. 

To estimate the photometric completeness, we put the 
1000 model quasars at each grid point in apparent magni- 
tude and redshift into Subaru Suprime-Cam FITS images 
as point sources, using the IRAF mkobjects task in the 
artdata package. As for 1000 model quasars, they were 
generated on a Monte Carlo method of drawing a value of 
alpha and EW. These point sources have apparent mag- 
nitudes calculated from their simulated spectra, in each 
image (g' , r', i', and z 1 ). We then tried to detect them 
in the z'-band image with SExtractor, and measure their 
colors in the double-image mode. Note that the mea- 
sured apparent magnitudes and colors of the simulated 
quasars are generally different from the magnitudes and 
colors before inserted into the Suprime-Cam images, due 
to the effects of photometric errors and neighboring fore- 
ground objects. Accordingly, some model quasars are not 
selected as photometric candidates of quasar with the cri- 
teria (1) - (5). To calculate the fraction of model quasars 
that are selected as photometric candidates in the above 
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Fig. 4. — Comparison between the empirical and simulated 
quasar colors. Solid black lines show the empirical mean colors 
of the SDSS DR7 quasars. Dashed black lines show 1 sigma dis- 
persions of the SDSS DR7 quasar colors. Red points show the 
simulated quasar colors. Red vertical lines show the standard de- 
viation of simulated quasar colors which do not include photometric 
errors. 



process, we estimate the photometric completeness at 
various magnitudes and redshifts (Fig. 5). The redshift 
range of the inferred completeness is moderately high at 
4-5 < z < 5.5. More specifically, the inferred complete- 
ness is ~ 90% for quasars with i' = 22.5, and ~ 80% for 
those with i' = 23.5 in that redshift range. The small 
dip at z ~ 5.2 — 5. 3 in the estimated completeness is due 
to the C IV emission that causes the i' — z' color to be 
red at that redshift range. 

5. QUASAR LUMINOSITY FUNCTION 

To calculate the upper limits of the quasar space 
density, we compute the effective comoving volume of 
the survey as: 

rz=co dy 

Ves(mi>) = (Kl C(m v ,z)—dz, (6) 

Jz=0 " z 

where dCl is the solid angle of the survey and C{m,i' , z) is 
the photometric completeness derived in Section 4. For 
comparison with other works, we convert the i'-band 
apparent magnitude to the absolute AB magnitude at 
1450 A (e.g., Richards et al. 2006; Croom et al. 2009; 
Glikman et al. 2010): 

Mi4 50 = m»' + 5 - b\ogd L {z) + 2.5(1 - a v ) 

log(l + z) + 2.5a„log (-^—) , (7) 



Fig. 5. — Simulated photometric completeness. Red, green, blue, 
purple, and cyan lines show the completeness for quasars with wiji 
= 22.0, 22.5, 23.0, 23.5, and 24.0, respectively. 

where di,(z), a v , and are the luminosity distance, 
spectral index of the quasar continuum (/„ oc v~ a "), 
and the effective wavelength of the i'-band (Ai'=7684A), 
respectively. We assumed the a v = 0.46 when we used 
the equation (7). As for the quasar candidates which 
did not perform the spectroscopic observations, we cal- 
culated the absolute magnitude at 1450A by assuming 
the effective redshift. Given the effective comoving vol- 
ume, the 1 a confidence upper limits on the space den- 
sity of type- 1 q uasars are calculated using statistics from 
IGehrelsi 111981 . 

The derived 1 a confidence upper limits on the space 
density of type-1 quasars are $ < 1.33 x 10~ 7 Mpc -3 
mag -1 for -24.52 < Mi 450 < -23.52 and $ < 2.88 x 
10" 7 Mpc~ 3 mag^ 1 for -23.52 < M li50 < -22.52. Note 
that there is another quasar candidate in the magnitude 
bin of —23.52 < M1450 < —22.52 which was not observed 
with FOCAS. We take into account the possibility that 
this candidate is a quasar when calculating the la confi- 
dence upper limit on the space density of type-1 quasars 
for -23.52 < M1450 < -22.52. The obtained 1 a confi- 
dence upper limits on the space density of type-1 quasars 
are plotted in Fig. 6. 

To compare our upper limits on the quasar space den- 
sity with the previous quasar surveys at similar red- 
shifts, we also plot the results of COMBO-17 (Wolf et 
al. 2003), SPSS (Ric hards et al. 2006), and GOODS 
(jFontanot et al. 20()3), in the redshift ranges of 4.2 < z 
< 4.8, 4.5 < z < 5.0, and 4.0 < z < 5.2 respectively, in 
Fig. 6. Note that the low-luminosity quasar sample of 
Fontanot et al. (2007) includes type-2 quasars while our 
sample and the sample of Wolf et al. (2003) do not in- 
clude type-2 quasars. To compare the result of Fontanot 
et al. (2007), we also calculated the quasar space den- 
sity when we included a type-2 quasar and the obtained 
quasar space density and its error are $ = 0.87~^q'j\ x 
10~ 7 Mpc" 3 mag" 1 for -23.52 < Mi 450 < -22.52 (see 
Fig. 6). This figure shows a marginal discrepancy be- 
tween the results of Fontanot et al. (2007) and of ours. 
However, the redshift difference between the two studies 
should be taken into account for such a comparison be- 
cause the quasar space density shows significant redshift 
evolution. 
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Fig. 6. — The z ~ 5 quasar luminosity functions. The red filed 
squares show our results (la confidence limits on the quasar space 
density) and the red open square shows the quasar space density 
when we include a type-2 quasar at z ~ 5.07. For clarifying the 
data plots in the figure, the open red square is slightly shifted to 
the left direction to avoid the overlap of the error bars. Triangles, 
squares, and circles denote the results reported by Fontanot et al. 
(2007), Richards et al. (2006), and Wolf et al. (2003), respectively. 

In Fig. 7, we plot the QLF reported by Fontanot et 
al. (2007) after correcting for the redshift difference (i.e., 
taking the redshift evolution in the QLF into account), 
adopting the model 13a in Fontanot et al. (2007). The 
model 13a assumes a pure density evolution of the QLF 
with an exponential form, that gives the minimum \ 2 
among the models examined in Fontanot et al. (2007). 
More specifically, in the model 13a, the bright-end slope 
is fixed to be 3.31 and there are three free parameters; 
are the faint-end slope, normalization, and redshift evo- 
lution parameter. Fig. 7 shows that our quasar space 
density at z ~ 5 are higher than the result of Fontanot 
et al. (2007) and thus our results are not contradictory 
to their result, once the redshift difference is corrected. 
Here it should be mentioned that Fontanot et al. (2007) 
adopted a different method in deriving the photomet- 
ric completeness from other surveys, that may introduce 
a possible systematic difference from other studies (see 
also Glikman et al. 2010). This effect is seen in Fig. 
6, where the inferred bright-end quasar space density is 
differen t between the r e sults by Fontanot et al. (2007) 
and by iRichards et al.l (|2006D even although the same 
SDSS quasar sample is used in the two studies. This 
suggests that the completeness adopted in Fontanot et 
al. (2007) may be underestimated, and accordingly the 
quasar space density is possibly overestimated by a fac- 
tor of ~ 2 — 3. In the case that we derive the faint end of 
the QLF at z ~ 5 by using the completeness which is cal- 
culated by Richards et al. (2006), the QLF of Fontanot 
et al. (2007) shifts toward lower space density in Fig. 7, 
i.e., well below our upper limits. 

In order to constrain the faint end of the QLF at 
z ~ 5 quantitatively, we search for parameter values that 
satisfy our result. Here we adopt the following double 
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Fig. 7. — Comparison of z ~ 5 quasar luminosity functions when 
we include a type-2 quasar at z ~ 5.07. The red squares show 
our results when we include a type-2 quasar at z ~ 5.07. The 
black dashed and solid lines are the QLF reported by Fontanot et 
al. (2007), before and after correcting for the redshift difference 
between their study and ours (assuming pure density evolution, 
see the main text for details). 

power-law function: 

$(M U50 ,z) = 

10 0.4(q + 1)(A^ 45 o-M* 450 ) + 1 o0.4(/3 + l)(A/ 145 o-Af- 450 ) >\°J 

where a, (3, $(M* 450 ), and M* 450 are the bright-end 
slope, the faint-end slope, the normalization of the lu- 
minosity function, and the characteristic absolute mag- 
nitude, respectively. Among the four parameters, the 
bright-end slope (a) is fixed to be a = —3.31 based on the 
SDSS results (see Fontanot et al. 2007). The parameter 
ranges which satisfy our results are shown in Fig. 8. Note 
that it is important to examine the redshift evolution of 
the faint-end slope and break magnitude, because such 
parameters give us useful constraints on the evolution- 
ary model of SMBHs and quasars (e.g., iHopkins et al.l 
2006, 2007). By taking into account of the results ob- 
tained in COSMOS for z ~ 5, the break magnitude in 
the QLF is brighter than M* 450 ~ —26 at z ~ 5. This 
is significantly brighter than the QLF break magnitude 
for z ~ 4 reported by Ikeda et al. (2011) and Glikman 
et al. (2011), as shown in Fig. 8. A possible explana- 
tion for this evolution is that the mass accretion in most 
quasars at z ~ 5 is higher than at z ~ 4 (although the 
quasar number density is lower at z ~ 5), which makes 
the characteristic magnitude brighter at z ~ 5 than at 
z~4. 

Here we discuss the evolution of the quasar space den- 
sity in the context of the AGN downsizing. The evo- 
lution of the quasar space density for different absolute 
magnitude bins provides important information to con- 
strain the evolution of SMBHs. Therefore we plot the 
quasar space density for different absolute magnitude 
bins as a function of redshift in Fig. 9. Although there 
are a number of l ow- luminosity quasar surveys at z ~ 3 
(jWolf et al.ll2TjQl iHunt et al.|[200l: iFontanot et aJ1l20"07t 
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Fig. 8. — Constraints on the faint-end slope and the break 
magnitude A^j* 450 . The shaded region shows the parameter space 
that is consistent with the inferred lcr confidence upper limits on 
the quasar space density at z ~ 5 derived in our study. A dot and 
triangle are the results inferred at z ~ 4 by Ikeda et al. (2011) and 
Glikman et al. (2011), respectively. Note that we estimate the er- 
rors of A4j" 45 q and the faint-end slope by applying a weighted least- 
squares method while Glikman et al. (2011) es timated them by 
apply ing the STY maximum-likelihood method (Efstathiou et al. 



iBongiorno et al.ll2007D . we plot only the results of the 
2dF-SDSS LRG and Quasar Survey (2SLAQ; Croom et 
al. 2009), the Spitzer Wi de-area Infrared E xtragalactic 
Legacy Survey (SWI RE: ISiana et all I2008TI. and SDSS 
([Richards et al.1 120061 ) . in order to avoid data with large 
statistical errors. While most studies at z < 3 suggest 
consistent results (i.e., the AGN downsizing), the situa- 
tion is rather controversial at z > 3. Recently, the new 
z ~ 4 QLF results of the NOAO Deep Wide-Field Survey 
(NDWFS) and the Deep Lens Survey (DLS) are reported 
by Glikman et al. (2011). Interestingly, the results of 
Glikman et al. (2011) suggest constant or higher num- 
ber densities of low-luminosity QSOs at z > 3, while our 
COSMOS result suggests a continuous decrease of these 
objects from z ~ 2 to z ~ 5. 

Our result is consistent with the downsizing AGN 
evolution suggested by previous quasar surveys at 
lower- z both in the optical and X-ray (e.g., Croom 
et al. 2009; Ueda et al. 2003; Hasinger et al. 2005). 
iWillott et all (|2010f ) reported the faint end of the QLF 
at z ~ 6 although there is a large error bar for the 
faintest magnitude bin because only one quasar was 
found. The z ~ 6 space density at -M1450 ~ —22 is lower 
than the upper limit of our z ~ 5 space density at the 
same magnitude and this result is also consistent with 
the AGN downsizing evolution. However, the results 
of Glikman et al. (2011) require a different picture at 
z > 3, being inconsistent with the downsizing scenario. 
It is not obvious what is causing this discrepancy. 
Masters et al. (submitted to the ApJ) stated that 
cosmic variance cannot be responsible for the observed 
discrepancy in space density of low-luminosity quasars 
between the COSMOS and the DLS NDWFS fields. If 
this discrepancy is due to the difference in the quasar 
selection criteria, then this discrepancy could be caused 
by the presence or absence of the point source selection 
on the HST images. However, both Ikeda et al. (2011) 
and Glikman et al. (2011) obtained the spectra of 
most quasar candidates to remove the contaminations. 




Redshift 



Fig. 9. — The redshift evolution of the quasar space den- 
sity. Red, green, blue, magenta, light blue, orange, black, 
purple lines are the space density of quasars with M1450 = 
-20,-21,-22,-23,-24,-25,-26, and -27, respectively. Filled 
squares, open circles, filled circle, and open triangles denote the 
results of different quasar surveys, as described in the upper-left 
side of the panel (see the main text for details). For the compar- 
ison, the results of COSMOS and previous surveys are connected 
by the dashed lines while the results of NDWFS k. DLS (Glikman 
et al. 2011) are connected by the dotted lines. 

Consequently, we cannot explain this discrepancy due to 
the difference of the quasar selection criteria. Therefore 
it remains important that we search for low-luminosity 
quasars at high redshift in other fields and derive the 
faint end of the QLF. We also plot the quasar space 
density measured in the GOODS fields (Fontanot et 
al. 2007) in Fig. 9, which is consistent with both 
results from Glikman et al. (2011) and COSMOS 
due to its large uncertainty. Further observations of 
low-luminosity quasars in a wider survey area are crucial 
to derive firm constraints on the scenarios of the quasar 
evolution, especially at z > 4. 

6. SUMMARY 

In order to examine the faint end of the QLF at z ~ 5, 
we select photometric candidates of quasars at z ~ 5 in 
the COSMOS field. The main results of this study are: 

• Although we discover the type-2 quasar at z ~ 
5.07, no typc-1 quasars at z ~ 5 are identified 
through the follow-up spectroscopic observation. 

• The upper limits on the type-1 quasar space density 
are $ < 1.33 x 10 -7 Mpc~ 3 mag^ 1 for -24.52 < 
M1450 < -23.52 and $ < 2.88 x 10" 7 Mpc~ 3 
mag" 1 for -23.52 < M 1A50 < -22.52. 

• The quasar space density and its error when we in- 
clude a type-2 quasar at z ~ 5.07 are $ = 0.87l ; 72 
x 10~ 7 Mpc~ 3 mag" 1 for -23.52 < M li50 < 
-22.52. 
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• The derived upper limits on the quasar space den- 
sity are consistent with the QLF inferred by the 
previous works at z <~ 5. 



The characteristic absolute magnitude of the QLF 
shows a significant redshift evolution between z ~ 4 
(M* 450 > -26) and z ~ 5 (M* 450 < -26). 



• A continuous decrease of the space density of low- 
luminosity (—24 < M1450 < —23) quasars is in- 
ferred, that is roughly consistent with the picture 



of the AGN downsizing evolution. 
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